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Abstract 

The diffraction of a plane wave by a transversely inhomogeneous isotropic nonmagnetic linearly 
polarized dielectric layer filled with a Kerr-type nonlinear medium is considered. The analytical 
and numerical solution techniques are developed. The diffraction problem is reduced to a singular 
boundary value problem for a semilinear second-order ordinary differential equation with a cubic 
nonlinearity and then to a cubic-nonlinear integral equation of the second kind and to a system of 
nonlinear operator equations of the second kind solved using iterations. Sufficient conditions of the 
unique solvability are obtained using the contraction principle. 



1 Introduction 



Scattering and propagation of electromagnetic waves in layered structures filled with nonlinear media have 
been a subject of intense studies since the 1970s. A goal of this work is to develop solution techniques for 
singular boundary value problems (BVPs) for the Maxwell equations arising in mathematical models of 
the wave diffraction in nonlinear media elaborated in [13] [16] , [27], [21], [22] that can be reduced to one- 
dimensional settings for the Helmholtz equation on the line [13] . [14] , The BVPs are formulated on infinite 
and semi-infinite intervals and with transmission-type conditions and conditions at infinity that contain 
the spectral parameter [27], [21], [22]. When the wave propagation in a cylindrical dielectric waveguide 
filled with a nonlinear medium is considered [21] , [22], the coefficient in the equation multiplying the 
nonlinear term differs from zero inside a finite interval (0, a) and the conditions are stated at the point a 
(continuity), at the origin (e.g. boundedness) , and at infinity (rate of decay). The corresponding singular 
semilinear BVPs are formulated for the differential operators L (A) u + a B (u; A) = 0, where Lu is a 
linear differential operator and B (u; A) is a nonlinear operator. An example with B (u) = u 3 associated 
with the study of the wave propagation in Kerr-type nonlinear fibers is considered in great detail in 
[21] . [22J. The method of solution employs reduction to nonlinear integral equations (IEs) [21], [22], 
[15] , [T6] . [25] constructed using Green's function of the linear differential operator Lu\ the eigenvalue 
problems are then replaced by the determination of characteristic numbers of integral operator-valued 
functions (OVFs) that are nonlinear both with respect to the solution and the spectral parameter. The 
latter problems are reduced to the functional dispersion equations, and their roots give the sought-for 
eigenvalues. The existence and distribution of roots on the complex plane are verified. The linearization 
is considered in [T5] . 

The reflection and transmission of electromagnetic waves at a nonlinear homogeneous, isotropic, non- 
magnetic dielectric layer situated between two linear homogeneous, semi-infinite media is of particular 
interest in linear optics |4J. In nonlinear optics, the Kerr-like nonlinear dielectric film has been the 
focus of a number of studies [8], [6], [7J, [ID], [12] • In [H], [6], [7], [10], the solutions of the nonlinear 
Helmholtz equations have been given in terms of various Jacobian elliptic functions. The explicit form 
of these functions depends on the associated parameter regimes. As shown in [12J, no classification of 
the solutions with respect to different parameter regimes is necessary, since the general solution can 
be presented in terms of Weierstrass' elliptic functions containing the complete parameter dependence. 
In [13], a simplified version of this result is given, generalizing the approach applied in linear optics. 
Namely, a general analytical solution of the Helmholtz equation is obtained describing the scattering 
of a plane, monochromatic, TE-polarized wave by a transversely homogeneous dielectric layer (with a 
constant permittivity) exhibiting a local Kerr-like nonlinearity. The layer is situated between two semi- 
infinite non-absorbing, non-magnetic, isotropic, and homogeneous media. The results derived contain the 
conditions for unbounded field intensities expressed in terms of the imaginary half-period of Weierstrass' 
elliptic function. The reflectivity R is calculated as a function of the layer thickness and the transmitted 
intensity. The critical values of R are determined. 

In [12] the approach set forth in [T3], [H] is applied to the analysis of the problems of the wave 
diffraction by layers filled both with linear and nonlinear dielectric media having constant and variable 
permittivities. The plane wave diffraction problem is reduced in [16J to a nonlinear Volterra IE and its 
solution is obtained as a limit of a certain function sequence. The sufficient conditions for the IE unique 
solvability are obtained by estimating the norms of the associated nonlinear operator. 

In this paper the approaches developed in [2TJ, [22] and [27], [32] [2D] are applied to the solution of 
singular semilinear BVPs arising in a mathematical model of the wave diffraction from a transversely 
inhomogeneous dielectric layer having a variable permittivity. The approach employs Fredholm-type IEs 
with complex-valued kernels derived on the basis of the method proposed in [27J and differs thus from 



the technique [13] . [16] . [25] based on the reduction to nonlinear Volterra IEs. On the other hand, the 
sufficient solvability conditions presented in this study are different from those reported in [16J; in fact, 
these conditions are obtained explicitly in terms of the problem parameters. Next, in this paper we apply 
the solution technique based on the analysis of cubic-nonlinear IEs to prove the unique solvability of the 
diffraction problem for a lossy weakly nonlinear layer with a complex-valued permittivity function. We 
note in this respect |25j where this problem is solved for a layer filled by linear and nonlinear lossy media 
using a general approach which enables one to evaluate the solutions in terms of uniformly convergent 
sequences of iterations of the Volterra IEs. 



2 Maxwell equations and wave propagation in nonlinear media 

2.1 General assumptions leading to the problem statement 

Nonlinear processes in electrodynamics and optics are described by the Maxwell equations 

VxE(f,t) = -I^l, Vx^(f, i) = ifl, 
V-Z>(r,t) = 0, V- B(r,t) = 0, 

Here E(f,t), H (f,t) D(r,t), and B (f,t) are the vectors of, respectively, electric and magnetic field 
intensities, electric displacement, and magnetic induction. This system is complemented by material 
equations 

D(f,t) = E(r,t) + 4irP{r,t), 
B (f, t) = H (f, t) + AtcM (f, t), 

— * —* 

where P(f,t) and M (r,t) are the vectors of, respectively, polarization and magnetic moment. 

The polarization vector P (r, t) = F E (r, t) , where F denotes a certain nonlinear operator, is 

generally nonlinear (with respect to the intensity) and nonlocal both in time and space In this work, we 
will limit the analysis, following [2J, to nonlinear media having spatially nonlocal response function. In 
this case the polarization vector can be expanded [3] in terms of the electric field components 

P t (f, t) = x^Ej + x%E 3 E k + x%E 3 E k E x + ... (3) 

Here Pi and Ei are the components of the polarization and electric vectors, respectively and coefficients 
X are lower terms of the expansion for nonlinear susceptibility. Thus, we ignore spatial dispersion pQ. 
Note however that this assumption does not limit our possibilities to consider effects in media with 
interfaces where linear and nonlinear susceptibility tensors may depend on spatial variables [2J. 

Below, we assume that the medium is nonmagnetic, M (r, t) = 0. Resolving equations ([1]) and (f2"l) 
with respect to H (r , t) we reduce them to one vector equation 

V 2 E (f, t) - V [V • E (f, t)\ - -^D^ (f, t) - ~T-^P {NL) (r, t) = 0, (4) 

where = E + 4nP^ = eE, P^ = x^E, and e = are the linear terms of the electric 

displacement and polarization vectors and permittivity tensor, respectively (here D\ L ^ = e^Ej, P^ = 

Xj] Ej, and e$ = 1 + 47ty,-^); p( NL ^ is the nonlinear part of the polarization vector (according to ©, 
pH - Xg&ft + W< + ..•); and*g>, ^ X g» , are the respeetive eomponents of the medium 
susceptibility tensors x^\ x^ ■, ■ 



Equation (@| is of general character and is used, together with material equations (j2j), in electrody- 
namics and optics; in every particular case, specific assumptions are made that enable one to simplify its 
form. Note, for example, that in the majority of important problems the longitudinal field components 



(along the z-axis) are negligible [2]. The second term in 



V • E (f, t) (where the inner product 



E- (Ve)/e 



can be written, using the conditionV • D = 0, in the formV • E = 

longitudinal and transverse field components and may be ignored in a number of cases 
Consider the case of stationary electromagnetic field 



contains both 



E (r, t) = Re exp (— iu t) E (f) 



H (r, t) = Re exp (— iu t) H (r) 



exp (— iu t) E {f) + exp (iu t) E* (r) 



exp (—iu t) H {f) + exp (iu t) H* (f) 



where E (r) and H (f) are the complex amplitudes of the electric, E (r,t), and magnetic, H (r,t), field 
intensity vectors, Re is the real part of the complex vector-function, and * denotes complex conjugation. 
Assuming that the medium is weakly nonlinear (when the so-called weakly-waveguide approximation 
holds), i.e. 

' « 4? , (5) 



SNL) 



where £^ L ^ = ^xfjii^kEl is governed by nonlinear terms in (jlj) (see [2], [26]), one can generalize 
weakly- waveguide approximation [21] and take into account the effect of nonlinear self-canalization [9]. 
In this case one can ignore the second term in (jlj), which is equivalent to ignoring the longitudinal field 
components, and vectors E and P will have only transverse components [2]. 

Consider the diffraction of a stationary electromagnetic wave [~ exp (—iu t)] by a weakly nonlinear 
object. Perform a transition to the frequency domain using the direct and inverse Fourier transforms 



.(3) 



E (r, u) 
(r,u) 

pi.NL) , p ~ s 



E{r,t) 
(f, t) 
. P^V (f,t) . 



e lut dt, 



E(r,t) 
(f, t) 

p(NL) ^ 



1 

2^ 



E (r, u) 
(f,u) 

p(NL) (f Cu] 



-i£jtj~. 



du. 



Applying formally the Fourier transform to equation (jlj) we obtain the following representation in 
the frequency domain 



V 2 E(f,u) 



,2 • 



+ ^LdM (f, u) + ^p(vl) {jf - ) = 

c c 



(6) 



A stationary [~ exp (—iu t)] electromagnetic wave propagating in a weakly nonlinear dielectric struc- 
ture gives rise to a field containing all frequency harmonics, see [T], [26]. Therefore, the quantities 
describing the electromagnetic field in the time domain subject to equation (jlj) can be represented as 
Fourier series 



E (r, t) — \ Yl E (f, nu) exp (—inut) , D( L > (r,t) — | (r, nu) exp (—inut) 

n=—oo n=— oo 

oo 

p(nl) 1) = i g pirn ^ nuj ^ exp (_ ina;£ ) . 



(7) 



Applying to (GO) the Fourier transform we obtain 



E (f, u>) 
(f,u 

p{NL) 



r,u 



I 



E(r,t) 
(f,t) 

p(NL) ( - £ ) 



n iCut 



dt 



I E 

-00 n=— oo 



P (r, ncu) 

£> (L) (r,W) 
p(7VL) ff* nuJ ^ j 



2tt 
2 



P (r, no;) 

D (L) (f, ncu) 
p(iVi) ^ nuo ^ j 



where 5 (s) = ^7= J exp (ist) dt is the Dirac delta-function. 



Substituting (jHJ) into (jHJ), we obtain an infinite equation system with respect to the sought-for Fourier 
amplitudes of the electromagnetic of the weakly nonlinear structure in the frequency domain, 

l 9 J 



V ■ E (r, nu) 



+ ^D {L) (f, ncu) + 



P (AfL) (f, ncu) 
0, ±1, ±2, 



0. 



For linear electrodynamic objects the equations in the system Q are independent. In a nonlinear 
structure, the presence of functions p( NL ^ (r, nu) makes them coupled since every harmonic depends 
on a series of E(f,nu). Indeed consider a three-component electromagnetic field E = (E x ,0,0), H = 
(0, H y , H z ). The fact that the field E = (E x ,0,0) has one component enables one to consider as a 
system of scalar equations with respect to E x . Take lower terms in the expansion in the vicinity of 
the zero value of the electric field intensity. Then the only nonzero component of the polarization vector 
P = (P,;, 0, 0) is determined by the third-order susceptibility tensor x^, which is characteristic for the 
Kerr- type medium. In the time domain, this component can be represented in the form (cf. ([3]) and ([7])): 



Px if, t) = i E p x (f, su) exp (-iust) = XxLxE x (f, t) E x (f, t) E x (f, t) 



P) 

s=— oo 

= | E Xxxxx {su; nu, mu, pu) E x (r, nu) E x (r, mu) E x (r, pu) e -^(n+ m + P )t_ (iq) 

j n, m, p, s = — oo 
[ n + m + p = s 

Applying to (fTUj) the Fourier transform with respect to time (jSj) we obtain an expansion in the 
frequency domain 



P x (f, su) = ~ XxLx (su ; nu, mu, pu) E x (r, nu) E x (r, mu) E x (f, pu) = 

j n,m,p— — oo 
| n + m + p = s 

oo 

= I E ZXxLx (su ; ju, -ju, su) \E X (r,ju)\ 2 E x (f, su) + 

3=0 

+ 4 E Xxxxx (su ; nu, mu, pu) E x (r, nu) E x (f, mu) E x (r, pu) . 

n, m, p = — oo 
n ^ —m, p = s 
m —p, n = s 
n ^ —p, m = s 
n + m + p = s 



The addends in the first sum of ( fill) are usually called the phase self-modulation (PSM) terms [2]. 
We obtained them taking into account the property of the Fourier coefficients E x (r,ju) = E* (r, —juj); 
factors 3 appear as a result of permutations {juj, —juj, suj} of three last parameters in the terms 

(3) i ' ' \ 

Xxixx [SUJ ; JUJ, —joj, SUJ). 

When particular nonlinear effects are considered, one can limit the analysis to finitely many equations 
of system (Q, leaving in the formulas for the polarization coefficients separate terms that characterize 
the physical problem in question. 

In this paper, we analyze electromagnetic fields scattered by a dielectric layer filled by a Kerr- 
type (weakly) nonlinear medium. We limit the analysis to such a level of intensities of the incident 
electromagnetic field affecting the structure when harmonic oscillations at combined frequencies may be 
neglected. In this case equations ([9]) and (fTTI) have the form 

V 2 E X (f, uj) - V [V • E x (f, uj)} + ^D X L) (f, uj) + ^P [ X NL) (f, u) = 0, 
P x (f, uj) = ~Xx%x {oj ; uj, -uj, uj) \E X (f, uj)\ 2 E x (r, uj) . 

In the next sections, we will use equations (fl2"]l to formulate a boundary value problem associated 
with the electromagnetic wave diffraction by a layer filled by a Kerr-type nonlinear medium and construct 
the methods of analytical solution to this problem. 



2.2 Statement of the problem of diffraction by a weakly nonlinear layer 

Denote by E (r) = E (r, uj) and H (f) = H (r , uj) the complex amplitudes of the stationary electromag- 
netic field; the time dependence is exp (—iut). Consider the problem of diffraction of a plane stationary 
electromagnetic wave E (r, t) = exp {—ioj t) E (r), H (r, t) = exp (— iuj t) H (f) by a nonmagnetic, M — 0, 
isotropic and linearly polarized E (f) = (E x (y, z) , 0, 0), H (r) = (o, H y = H g = -^^) 

(E-polarization), transversely inhomogeneous, e^(z) = e x L x \z), dielectric layer with a weak Kerr-type 
nonlinearity §1) P ( X NL) = (3/4)x£Lx \E X \ 2 E x (where P^ NL > = (P ( X NL \ 0, 0); this problem is stated in 
[IS], PJ, P], and [25 . Using (□), (EJ, and the results from [2] we obtain V • E = -E (Ve)/e from the 
equation V-D = 0; therefore, the second term is absent, both in (j3J) written in the time domain and in 
©, &, and (H2D, V (V-E\ =0. 

According to [18J and the results of the previous section, one can show that the total field E x (y, z) = 
E l x nc (y, z) + E^ cat (y, z) of diffraction of the plane wave 

E™ c (y,z) =a inc exp{i [0 y - Y (z - 2 tt 5)} } , z > 2 tt <5, 

by a weakly nonlinear dielectric layer is the solution to the equation (see (j4])): 

V 2 ■ E + ^) (z)-E+ *™LpW) = [v 2 + k 2 b (z, a, \E X \ 2 )} E x (y, z) = (13) 
satisfying the following generalized boundary conditions: 

continuity of E tg and H tg on the boundary of the layer having the permittivity e (z, a, \E X \ 2 ) , 
the spatial quasi-homogeneity condition [27] with respect to y 



E x (y, z) = U (z) exp (i<j> y) 



(14) 



and the radiation condition for the scattered field 



n scat ( \ 

E s x cat (y,z) = l bscat | exp(z {(f) y ± T {z T 2 ti 5))) , z < ±2tt5. (15) 

Here we use the following notations: {x, y, z, t} are dimensionless spatial-temporal coordinates intro- 
duced so that the layer thickness is Air 5; the time dependence is exp (— % to t); uo = ncis the dimensionless 
circular frequency; k = uo/c = 2ix/\ is the dimensionless frequency parameter such that h/X = 2k S, 
where A is the free-space wavelength; c = (e / i o) 1//2 is the dimensionless quantity equal to the speed 
of light in the medium containing the layer (Imc = 0); Eq and /io are the material parameters of the 
medium; E tg and H tg are the tangential components of the total E and H fields; V 2 = d 2 /dy 2 + d 2 /d z 2 ; 



e (z, a, \E X \ 2 ) = e (z, a, \U (z)\ 2 ) = j 



1 , \z\> 2n5 
eM(z) + a\U (z)\ 2 , \z\<2n5 



where (z) is piecewise continuously diflerentiable with respect to z\ a = 37rXxxxx] T — (k 2 — <p 2 ) x l 2 
is the transverse propagation constant (transverse wavenumber); <f> = k • sin(</?) is the longitudinal 
propagation constant (longitudinal wavenumber); and tp is the angle of incidence of the plane wave, 
\ip\ < 7r/2. Quantities x' ,y' , z' ,t' ,u' are reconstructed from the dimensionless values by the formulas 
(x 1 , y', z') = (x,y,z) • h/AitS, t' = t ■ h/An 5, and uj' = uin 5/ h. 
We look for the solution to problem (jlcfl) - (115p in the form 



E x (y, z) = U (z) exp (i <j> y) = 

a inc Qxp |- ^ y _y(z- 2ti5)}} + a scat exp {i [<py + V (z - 2tt5)}} , z > 2tt5, , . 

U scat (z)exp(i(f)y) , \z\<2n5, ^ > 

h scat exp ^ ^ y _y( z + 2tt5)} } , z < -2tt5, 

assuming the continuity on the permittivity break lines z = 2%8 and z = —2tt5, so that U (—2ir5) = b scat 
and U (2n5) = a inc + a 



3 Integral equation of the nonlinear problem 

3.1 Reducing to an integral equation 

We solve problem f|T3|) - f|T5l) in the whole space Q = {q = {y, z} : — oo < y, z < oo} by reducing it to a 
one-dimensional IE along the layer height z G [— 2tt5, 2ix5\ with respect to the scattered field component 
U (z) = JJ scat (z) introduced in f|T6|) . To this end, make use of canonical Green's function Go of problem 
(H3D-(H5D (for e = l) defined in the strip Q { y j0o} = {q = {y, z] : -Y <y <Y , \z\ < oo; Y > 0} C Q 
by the expression [TS], [T7j, [T9] 



Go (q , 9o) = Ty exp {i [(f) {y - y ) + T \z - z \] } /V = 



exp(±z4>y)§ J tf 



(i) 



K 



\j{y- vo) 2 + ( z ~ z oY 



exp (Ti<fry) dy . 



(17) 



The nonlinear IE with respect to U (z) introduced in (TT6l) is obtained using a classical approach set 
forth in [23J. Denote by V (q) = E x (q = {y, z}) = U (z) exp (i<fiy) the total diffraction field (see f|T6|) ), 
where U (z) is the solution of problem f|T3|) - f|T5|) . and write equation f|T3|) as 



(V 2 + k 2 ) V(q)=[l-e(q,a, \V (q)\ 2 )} k 2 V (q) 



The field of the incident plane wave Vo (q) = V mc (q) = a mc exp {i [<fiy — T (z — 2tt5)] } satisfies in the 
whole space Q the homogeneous Helmholtz equation 

(V 2 + k 2 ) V (q) = 0. (19) 

At z > 2tt5, Vq (q) is the incident field of the incoming plane wave irradiating the layer, while at 
z < 2n5 Vo (q) is the outgoing plane wave that satisfies the radiation condition at infinity (because in 
the representation for Vo (q) the transverse wavenumber T > 0). 

Subtracting from (IT51) equation (TT9l) we obtain 

(V 2 + k 2 ) [V(q)-V (q))=[l-e(q,a, \V (q)\ 2 )] k 2 V (q) . (20) 

Here V (q) — V (q) satisfies the radiation condition (fl5j) in the whole space. In fact, at z > 2tt5 the 
difference V (q) — Vo (q) = V scat (q) is the reflected field and at z — ► — oo both V (q) and Vo (q) satisfies 
the radiation condition. 

Using (120]) and the equation for the canonical Green's function Go 

(V 2 + K 2 ) G (q,q ) = -5(q,q ) (21) 
(where 6 (q, go) is the Dirac delta- function) it is easy to show that 

(V - Vo) V 2 G - G V 2 (V -V ) = -(V- V ) 6 (g, q ) - G [l - e (q, a, |V| 2 )] k 2 V. (22) 

Let Q{y,z} = {q = {y, z} '■ —Y < y <Y, —Z < z < Z; Y > 0, Z > 2n5} denote a rectangular do- 
main in space Q. Divide this domain into rectangles such that in each of them the permittivity 
e (q, a, \V\ 2 ) is continuously differentiable with respect to y and z. On common parts of the boundaries 
of rectangles V (q) and dV (q)/dn (where n denotes the outer normal) are continuous due to the conti- 
nuity of the tangential components E tg and H tg . Therefore, in the whole domain Q{y, zy, the sought-for 
twice continuously differentiable function V (q) preserves this property up to the boundary 0Q{y, z}', i-e., 
V (q) E C 2 (Q {Y ,z } ) n C 1 (Q{y,z } ) (here Q {Y , Z} = Q {Y ,z } U dQ {Y , z} ). 

Applying in Q{y,z} Green's formula 




[(V - Vo) V 2 G - G V 2 (V - Vo)] dq = j 



Q{Y,Z} 

and taking into account 



9Q{y,z) 



(V - V ) 



dG 
dn 



G, 



d(V~V ) 
dn 



dq Q 



, we obtain 



V{q) = -K 2 J f G (q,q )[l-e(q ,a, \V (q )\ 2 )} V (q ) dq + V (q) - 

Q{Y,Z} 

I {[V(q )-V (q )]^t^-G (q,q ) 



9Q{y,z} 



d[V(q )-V (q )] 
dn 



dq Q 



(23) 



When parameter Z — >• oo, the integrals in the lower, [(—Z, —Y) , (-Z, Y)], and upper [(Z, Y) , (Z, —Y)] 
parts of the boundary 8Q{y,z} that enter curvilinear integral ( 1231) tend to zero. This statement follows 
from asymptotic properties of Green's function ( flOl) and the fact that V scat = V — Vo satisfies radia- 
tion condition ( Tl5j) . The integrals along [(— Z, Y) , (Z, Y)] and [(Z, —Y) , (— Z, —Y)] cancel each other. 
Therefore, setting in (f23|) Z — > oo and deleting the curvilinear integral along the boundary 8Q{y,z^oo}, 
we obtain an integral representation for the total field of diffraction in the band Q{y,oo} 



V(q) = - 




G (q, q )[l-£ (g , a, \V (g )| 2 ) ] V (q ) dq + V (q) , g G Q { y i00 }. 



The integrand in the remaining double integral is a finite function with respect to z, i.e. 

l-e(q,a,\V(q)\ 2 )=0, \z\ > 2tt<5, 

because the permittivity of the medium enveloping the layer is assumed to be equal unity, so that one 
can limit the integration to the domain occupied by the dielectric 



V(q) = -K 2 j J Go(q,qo)[l-e(q ,a,\V(q )\ 2 )]V(q )dq + Vo(q), q G Q{y,oc}- 



Q{Y,Z=2wS} 



Performing a transfer to the limit Y — > oo (which can be justified by the facts that, according to 
(j!4p and ffTTl) . the integrand is asymptotically equivalent to O (Y^ 1 ) and parameter Ym&y be chosen 
arbitrarily) we obtain an integral representation for the total field of diffraction in the whole space Q 



V(q) = -K 2 J J G (q,qo) [l - e (go, a, \V (q )\ 2 )} V (q ) dq + V (q) , q e Q . 



(24) 



Here Qs = Q{oo,z=2-kS} = {<? = {y, z} : — oo < y < +00 , \z\ < 2n5} is the band occupied by the non- 
linear dielectric layer. 

We can also obtain using an iteration scheme based on the approach developed in [IB] , [H] . Let 
us give a short description of this method. In space Q a function sequence V n (y,z) is constructed such 
that every function of this sequence, beginning from n = 1, satisfies conditions (TUT) and (fT5l) . and the 
limiting function V = E x (y, z) = lim V n is a solution to (|T3|) - (|T5|) ; namely, 

(V 2 + k 2 ) V = 0, (V 2 + /t 2 ) V 1 =[l-e(z, a, \V Q \ 2 )] k 2 V + V , . . . , . , 

(V 2 + k 2 ) V n+1 = [l-e(z,a, |K| 2 )] K 2 V n +V 0} ... 1 ' 

Equations ( J25l) are formally equivalent to the following 



V (q) = V ins (q) , 

Vi (q) = -JJG q (q, q ) [l - e (q , a, \V (q )\ 2 )] k 2 V q (q ) dq + V (q) , . . . , 

Qs (26) 

Vn+i (q) = - II G (q, q Q )[l-e (q , a, \V n (q )\ 2 )] K 2 V n (q ) dq + V (q) . . . , q £ Q. 

Qs 

Performing in (1261) a transfer to the limit n — * 00 we obtain the integral representation (j24p for the 
total field of diffraction in Q. 

For q G Qs, representation ff24l) is transformed to a nonlinear IE with respect to the sought for 
scattered field V (q) = V scat (q) , q G Qs, see f|T6|) . Substituting into equation fT28|) formula ffTTj) for 
canonical Green's function and the expression for the permittivity e (qo, a, \V(qo)\ 2 ) we obtain an 



equation 

U(z) e i4,y = 

{2n8 Y \ 
J e^V 1 ^ [1 - (e<*> («,) + a |f/(^o)| 2 )] U(z )dy dz 1 (27) 

+ t/ inc (2) e'*" 

with respect to C/(z) = £/ scai (z) , \z\ < 27r<5 , that enters the expression for the field V (q) = 
(g = {$/, 2;}) = U (z) exp quasi-homogeneous along the layer. Integrating in domain Qs with 
respect to yo, we obtain a nonlinear IE of the second kind with respect to an unknown function 
U(z) E L 2 ([-2ir5, 2tx5}): 

2tt<5 
2 r 

U{z) + ^ J ex P (ir\z-z \)[l-(e^(z )+a\U(z )\ 2 )]U(z ) dz (28) 

-2tt5 

= U inc (z), \z\< 2n5, 

where U mc (z) = a mc exp [— iV (z — 2nS)]. 

The existence and uniqueness of the solution to IE (128]) for the linear problem with a = are proved 
in [18], [28]. In the general case a nonlinear IE of type (1281) may or may not have (the unique) solution. 
Its solvability is governed by the properties of the kernel and the right-hand side (incident field U mc (z)) 
and value of the nonlinearity parameter. 

Note that from the method of obtaining IE (T28]) it follows that the solution to this IE may be used for 
the integral representation in E x (y, z) = U (z) exp (icjyy) of the sought for solution to problem (fl3~|) - (|15p 
for the points with the coordinates \z\ > 2tt5 outside the nonlinear layer. Indeed, finding the solution to 
IE ( 12"B"|) and substituting it under the integral sign in (|28[) . we obtain an explicit expression for U (z) in 
the domain \z\ > 2tt5. 

The equivalence of IE ([28"]) to problem (|T3]) - (fT6]) is proved in Appendix. 

3.2 Sufficient condition of the existence of solution to nonlinear IE 

Assume that the permittivity function (zq) is positive, bounded, and continuous in the interval 
7 = [— 2it5, 2tt8], so that (z ) E C (7), where C (7) denotes the space of continuous functions in the 
closed interval 7 with the norm ||/|| = H/IU^ = max \f (z)\. Assume also that 

1 < e {L) (z) < E , z €7, E>1. (29) 
Write IE ( 128]) in the operator form 

U + AU - a F (U) = f , (30) 

where 

/ ( z ) = U inc (z) = a exp [-i k cos (ip) (z - 2ir5)} , a = a inc > , 



2tt8 

AU= J k(z- z ) [1 - (zq)] U (z ) dz 

-2ttS 



(31) 



is a linear integral operator with the continuous kernel 



% tv /IT 7T \ 

k(t) = s exp[2iK cos (tp) \t\) , s = - — \2 <V< 2J' 

and 



2tt<5 



F(U)= J k(z-z )\U(z )\ 2 U(z ) dz 

-2tv8 

is a cubic-nonlinear integral operator. 

A linear integral operator B : C (7) — > C (7) defined by 

2tt5 

£>w = J k (z — z )u (z ) dzo 

-2tt5 

is bounded, so that 

2tt(5 



\B\\< I max\k(z - z )\dz = 4n5- — = — = q . (32) 

2 cos (</?) cos (93) 



-2tt<5 



Integral operator (1321) is bounded and continuous in C (7) and its norm can be estimated as 

2w8 

\k(z-zo)\ |1-£ (L) (z )\dz 



< max 

267 



-27n5 



£-1 __= E-l)^ (33) 
cos (99) 



The nonlinear operator Q (U) = \U\ U is bounded and continuous in C (7) and F (U) is therefore 
bounded and continuous in C (7) as a superposition of .£> and Q. Hence the nonlinear operator T (U) = 
—AU + a F (U) + / is completely continuous on each bounded subset fl G C (7). 
Set 

n = s p = {uec( 1 )-. \\u\\< P } 

to be a ball in C (7), assume that Z7 G 5* p , and estimate the C (7)-norm of T: 

||T(f/)|| < \\A\\ \\U\\+a \\B\\ \\Uf + 11/11 < 



< p \\A\\ +a \\B\\ p 3 + a < K (p) , U G S. 



p ■ 



K (p) — (E — 1) g P + a <?oP + a 



Write equation fl30l) as U = T (U). One can apply to the operator T ([/) : 5* p — > S p the Banach fixed- 
point theorem [11] if K (jp) < p . In order to determine the corresponding range of values of parameter 
p > 0, solve the inequality K (p) < p , which yields a cubic inequality with respect to p 

P K (p)=D p 3 -D lP + a<0, D = aq , D l = 1 - (E - 1) q . (34) 

The necessary condition for ( 1341) to have a positive solution is .Di > 0, which yields (E — 1) qo < 1, 
or, according to (|32|) and assumption (|29|) concerning the properties of the permittivity function (z), 



2tt5k < cos (y?) <^ max [e {L) (z)] -1} (--< ip < -) . (35) 



Subject to the condition (l3~5j) . cubic polynomial Pk (p) in ( 1341) has two positive zeros p\ and P2 if 

< a < max(— Dop 3 + D\p) (with Z)i > 0), which holds if the local minimum of Pk{p) is negative, 

p>0 

Pk (Pext) < 0, at the point p ext = J > where P' K (p ext ) = 0. The corresponding condition for a can 
be written as a < |Diy / ^^-, or 

2 [! — (£ — 1) q f 2 2n5K 

3V3 J% cos (<p) 



Condition (I36p holds for arbitrary set of the problem parameters a, k, ip, 5, and E satisfying (1351) if 
the nonlinearity parameter a is sufficiently small because a enters only the left-hand side of inequality 
(1361) . The inequality K (p) < p holds forp G (pi, P 2), Pi > P2 > 0; for example, at 



(E - 1) g 



3 a go 



V = Pext = 

Theorem 1 Assume that 

(i) the permittivity function e^ L > (z) is positive, bounded, and continuous in the closed interval 7 = 
-27r<5, 2vr<5] and E = max (z)] > 1; 

(ii) parameters a, k, 5 and the nonlinearity parameter a are positive, \(p\ < n/2, and all the parameters 
a, k, 5, ip, a, and E satisfy §E2j) , (dSP and l[3b}) . namely 



j-, 1 / j-, 1 \ , 4 1 [1- (£-l)g ] 3 2ttSk 
E>1, (£-l)g <l, a<a = —— v - * , g - 



27 a 2 go cos (cp) 

Then the operator T(U) = -AU + a F(U) + f, T (U) : S p -»• S p defined by $® and (TUP is a 
contraction in the space C (7) if 

t = q {E-l + 3ap 2 ) < 1, p G ( Pl ,p 2 ) (37) 

where p\ and p 2 are positive zeros of the polynomial Pr{p) defined by ( Tff^j ) and a is sufficiently small, 
satisfying 

< a < min{a ,ai} , ol x = - {g [ 1 - (E - 1) g ] } _1 . (38) 

Proof. Use definition fl30|) . estimates fl32l) and (1331) . and inequality | |zi| — j^l | < |^i — z 2 | (where 
zi, z 2 are complex numbers), assume that U, V G 5* p , and estimate the C (7)-norm of the difference 
r(C7)-T(^): 



\\T (U) - T (V)\\ < \\AV - AU\\ + a \\F (U) — F (V)\\ < 

< (E -l)q \\U -V\\+aq \\U\U\ 2 -V\V\ 2 \\ < g (E - 1 + 3ap 2 ) \\U-V\\ . 

Thus, inequality (!37|) provides that operator T ([/) : — > is a contraction if a is sufficiently 
small; namely, satisfies (|36|) and (I38I) . Note that (|38|) follows from (I36I) . the condition 3apg Xi < 1, and 
inequality < p\ < p ext , where p ext G ij>i,p 2 ) is the point of a negative local minimum of the cubic 
polynomial Pk ip) (1511) satisfying P' K {p e xt) = and mm [Pk if)} = Pk ip e xt) < and pi, p 2 are positive 

zeros of Pk (p) ■ 

Summarizing the results verified above we conclude that the following statement is valid. 



Theorem 2 Assume that the permittivity function e^ L > (z), parameters a, k,5, the nonlinearity parameter 
a, and quantity E satisfy conditions (i) and (ii) from Theorem 1 and conditions (37) and (QSP- Then 



the operator T (U) = —AU + a F (U) + f defined by fifty and fiJI\) is a contraction in the space C (7) 
and IE fifty has the unique solution U* (z) continuous in the closed interval [— 2tiS, 2n8\. U* (z) is a 
limit with respect to the C (7) -norm of the function sequence U n (z) (the fixed point of operator T (U)) 
determined according to 

U n+1 = T{U n ), 7i = 0,l,2,..., U ES P = {U EC^): \\U\\<p}, p 6 (pi,p 2 ), (39) 
where p\ and p 2 are positive zeros of the polynomial Pk (p) defined by fiffy . 

The rate of convergence of the fixed-point iterations (1391) can be estimated using the quantity t < 1 
defined in (13] 



\\U n -U*\\ = \\T(U n ^)-T(U*)\\ <t \\U n ^-U*\\ <...<t^ L \\T(U )-U*\\, 71 = 2,3,... . 

3.3 Complex-valued permittivity function (diffraction by a lossy nonlinear 
layer) 

The method and results can be extended to the case when the permittivity (z) is an arbitrary 
complex-valued function of the real argument z continuous and bounded on the line. To this end denote 



: {L) (z) - 1 



g [z) = e 1 (z) exp [ie 2 (z)\ = g x (z) + % g 2 {z) 



(40) 



where, according to physical assumptions of the model, the real and imaginary parts of the permittivity 
function, g\ (z) and g 2 (z), are positive, continuous, and bounded on the line satisfying g\ [z) > 1 and 
g\ (z) » g 2 (z), so that the modulus E\ (z) and argument e 2 (z) of the permittivity are also positive 
functions continuous and bounded on the line with < e 2 (z) < tt/2. 
Make use of (HUj) and represent integral operator (|3T|) as 



2ti-<5 



AiU 



k x (z, z ) £1 (zq) U (z ) dz , kx (z,z ) = -s exp{i [2ncos(ip) \z - z \ e 2 (z )}} • (41) 



-2tt5 



Assuming, similar to (|29~!) and taking into account (|40j) and the conditions for the permittivity func- 
tion, that 



< ex (z) < Ex , z e 7, (42) 

(that is, < (z) \ < Ex, z G 7) we can estimate, as in ( J3T|) . the norm of the integral operator (HI]), 
which is bounded and continuous in C (7), as 



|| < max 

267 



2n8 



\kx (z,z ) \ \sx (z )\ dz 



-2tt5 



2ty5k 
Ei t— r = Ex go- 
cos (ip) 



(43) 



3]) yields an estimate for the norm of the nonlinear operator Tx (U) = —Ax U + a F (U) + / 



\\Ti(U)\\ < K x {p) , UeS p , K 1 (p)=E 1 q p + aq p 3 + a. (44) 

Thus one can easily check that the following statements are valid which are extensions of Theorems 
1 and 2 to the case of a complex-valued permittivity function. 

Theorem 3 Assume that 

(i) the permittivity (z) is a complex-valued function given by ( f^ffi ), where g\ (z) > and are contin- 
uous and bounded on the line so that the modulus £\{z) and argument £2(2) of the function (z) - 1 
are also nonnegative functions continuous and bounded on the line with < £2 (z) < tt/2 and £\ (z) > 1, 
and Ex = max [si (z)] > in the closed interval 7 = [— 2n5, 2uS\; 

(ii) parameters a, k, 5 and the nonlinearity parameter a are positive, \tp\ < tt/2, and all the parameters 
a, k, 5, ip, a, and E\ satisfy the conditions similar to (f^j. |23J) and l[3(j\) . namely, 

m 4 1 (1 -E iqo ) 3 2tt5k 

0<E iqo <l, a<a^ = —-± q = > 0. 45 

27 a 2 q cos {(p) 



Then the operator Ti (U) = —A\U + a F (U) + f, T x (U) : S p — > S p defined using ( [^7| j is a 
contraction in the spaceC (7) if 

h = q (E 1 + 3ap 2 ) < 1 , p G (p?^) , (46) 
where p^ and p^ are positive zeros of the polynomial 

pW(p) = D p 3 -D^p + a, D = aq , = 1 - E 1 q (47) 

and a is sufficiently small, satisfying 

0<a<min {a®, a®}, aj 3 = ^ [go (1 - ^o)] -1 ■ (48) 

Theorem 4 Assume that the permittivity function £^ L > (z) specified by $JW , parameters a, n, 5, the 
nonlinearity parameter a, and quantity E 1 = max [e\ {z)\ > 0, 7 = [— 2nS, 2ttS] (e\ (z) = |e^ L ^ (z) — 1 ) 

satisfy conditions (i) and (ii) from Theorem 3 and conditions $45\ ), [4°V > an d ( fT^P - Then the operator 
T\ {U) = —A\U + a F (U) + f defined using is a contraction in the space C (7) and IE $28\) has the 
unique solution U* (z) continuous in the closed interval [— 2tt5, 2tt5]. U* (z) is a limit with respect to the 
C {^j)-norm of the function sequence U n (z) (the fixed point of operator T\ (U)) determined according to 

U n+1 = Ti (U n ) , n = 0, 1, 2, U G S p = { U G C (7) : \\U\\ < p} , p G (p? , p^) , (49) 
where p^p and p^ are positive zeros of the polynomial Pjp (p) defined by ( f^7| ). 



The rate of convergence of the fixed-point iterations (1491) can be estimated using the quantity ti < 1 
defined in (1461): 



|| U n -U*\\< fr 1 II T x (U ) -U*\\, n = 2, 3, . . . , 
U e S p * = { U G C (7) : \\U\\ < p*} , p* G (p?,^ 



The existence of the unique solution to IE (|28l) subject to the sufficient conditions specified in for- 
mulations of Theorems 2 and 4 (corresponding to the cases of, respectively, real- and complex-valued 
permittivity function of the nonlinear layer) and the equivalence of IE (1281) to problem ( 1X31) ( TT6|) proved 
in Appendix enables us to prove the following statement which constitutes the main result of this study 

Theorem 5 Assume that the permittivity function (z) is (a) real-valued and quantity Eand param- 
eters a = a inc , k, 5, if, and a satisfy conditions (i) and (ii) from Theorem 1, UTV), and or (b) 



complex-valued (given by KJW) and quantity E\ = max (z) — 1 > and parameters a = a inc , n, 

zS7 

5, ip, and a satisfy conditions (i) and (ii) from Theorem 3 and conditions fij5\ ), jjj^ )> an d ( fT^P - Then 
problem (T3})- (Tb)) has the unique solution U* (z) continuous in the closed interval [— 2n5, 2n5] which 
can be determined as a limit with respect to the C ('j)-norm of the function sequence U n (z) determined, 
respectively, according to (a) |22)) or (b) ( fT^P - 

3.4 First iterations as trigonometric polynomials 

In view of the fact that at (z) = 1 and a = nonlinear IE ffTBl) has a formal solution 



U (z) = U mc (z) = a exp (— i bz) , a = exp (ibd) , d = 2n5 , b = k cos (<p) , (50) 

it is reasonable to choose the zero iteration in ( f49i) Uq (z) = a exp (— ibz) in the form ( 1501) . The linear 
integral operators (1311 and (HTl) can be represented as 

d 

AU = A[rj\U = J k(z-zo) J] (zb) U (z ) dz , 
k(t)=8 exp{2ib\t\), s = T ^, (-|<^<f). 

Obviously, they are linear with respect to the (continuous complex- valued) weight function rj (zo): 

A [hirji + h 2 rj 2 ] U = hiA [771] U + h 2 A [rj 2 ] U , hi, h 2 = const. 



Lemma 1 

where 



A [r/ 0) ] U = Hi exp (-z 2bz) + H 2 exp (z 2bz) + H 3 exp [iz{q- b)} 

iT osq 



Hj = H H 3 , (j = 1, 2, 3) , H { 







(b + q) (Sb-q)' 



Hi = (3b-q)exp[id{b + q)} , H 2 = {b + q) exp [i d (36 - q)} , H 3 = -46, (52) 
q 7^ —6, q 7^ 36, and 

U (z) = a exp (—ibz) , r/ -* (2) = Texp (iqz) , a, T — const. (53) 

At g = 6 ; 

A ^ (0) ] C/ = [exp (i 2bd) cos (z26^) - 1] . 



Proof. Proof of Lemma 1 reduces to tedious algebra and integration. 

We see that it is possible to determine explicitly the image A M )] Uq of a simple trigonometric 
polynomial Uq = aexp (— ibz)a.nd to show that this image is also a trigonometric polynomial: 

3 

A [V - 1 ] [a exp (— % bz)] = Hj exp (i Cjz) , c\ = —2b , c<i = 2b , c 3 = q — b. 

The linearity of A [rj] U with respect to the weight function rj (z) and U yields 
Lemma 2 Let 

(z) = rj exp (i qjz) , > 1. 
3=1 

Then the image A [?/ )] U of a trigonometric polynomial U (z) = ^ hjexp (ibjz) is also a trigono- 
metric polynomial: 

n a 

A[r)W] [/(*) = ^ P3 exp (ie^), 
3=1 

where the coefficients Pj and the number of terms can be determined explicitly. 

Similar statements are valid for the nonlinear operator F (U) defined in (I3"0l) ; namely, the image 
F(Uq) of a trigonometric polynomial Uq = aexp (— ibz) is also a trigonometric polynomial that can be 
determined explicitly: 



Lemma 3 

where 



F (Uq) = fi exp (— i 2bz) + $2 exp (i 2bz) + f 3 exp (— i bz) 



f — f f f„-_19Q\ f _ ia 3 exp(ibd)s 

Jj — JoJj , [J — 3 3) ) Jo — 3^ ' 

/1 = 3 exp (idb) , f 2 = ~ exp (idb) , f 3 = -2 , 

[To (-2) = 5 exp (— z te) , a = a exp (2 6d) , a = const . 

Lemmas 1-3 enable one to evaluate explicitly the first iteration 

U\ = T (Uq) = -AUq + aF (U ) + U = 
3 2 
= — Hj exp (i cjz) + a^] fj exp (i cjz) + af 3 exp (— i bz) = 
3=1 ' 3 1 

4 

= Sj exp (i Cjz) , S'j = -if, + afj (j = 1, 2) , S 3 = af 3 , S± = -H 3 , 
3=1 

c\ = —2b , c 2 = 26 , c 3 = —b , C4 = q — b . 

We conclude that according to ( 13T|) and ( 15T1) if the permittivity function e^ L ^ (2) is a trigonomet- 
ric polynomial then the first iteration U\ specified by (1491) is also a trigonometric polynomial whose 
coefficients can be determined explicitly. 



3.5 Sufficient condition of the existence of solution to nonlinear IE: reduc- 
ing to a functional equation system 

Here we present the proof of an alternative sufficient condition for the existence of a solution to nonlinear 
IE ( 1281) which is similar to the solvability conditions of the type ( 1371) . The approach developed in this 
section enables one to create a rather efficient method of the numerical solution of the IE. To this end, 
reduce ( 128]) to a nonlinear functional equation system, considering the system of two IEs in the domain 

\z\ < 2n5 : 

U n+X (z) + l 4 Y exp (iT \z - *o|) [1 - (eW (z ) + a \U n (z )\ 2 )] U n (z ) dz = U mc (z) , 

-?* (54) 
m n {z) + l 4 j exp (iT \z - z \) [l - (e (L) (zq) + a \U n (z )\ 2 )] # n (z ) dz = U mc (z) . 

-2tt<5 

The first equation of system (1541) is an iteration scheme of solution to nonlinear equation (f28l) 
(cf. (1261)). The second is a linear IE with respect to (z) for the given U n (z ) . If \l/ n (z) is not 
an eigenfunction of the problem of diffraction by the layer with the permittivity e (z, a, \U n (z)| 2 ) = 
(z) + a \U n (z)\ 2 , then the second equation is uniquely solvable [15, 26] and its solution can be 
represented in the form 

tf ft (z)=*(z,a, \U n (z)\ 2 )U mc (z). (55) 

where \1/ (z, a, \U n [z)\ 2 ) is the solution to the linear IE at U mc (z) = 1 such that (z, a, 

The analysis of the convergence criterion for the sequence U n (z) , \l/ n (z) specified by system ( 15"4I) 

enables one to obtain a sufficient condition for the existence of solution to nonlinear IE (1281) . 

Kernels of IEs (154]) are identical, which makes it possible to calculate and estimate the L2-norm of 

the difference between U n (z) and \l/ n (z) 



p\U„ 



z) , ty n (z) 



n+l \Z) , 



2tt<5 



/ \U n+1 (z)-V n (z)\*dz 



-2tt<5 



1/2 




2ir8 



J exp (zT \z - z \) [1 - (e (L) (z ) + a\U n (z )\ 2 )] [U n (z ) - * n (z )} dz 

-2ir5 



1/2 



dz 



2tt& 



J [1- (e^(z ) + a\U n (z )\ 2 )} [U n (z )-V n (z )]dz 

~2tt8 



1/2 



< 



<*1 
— 2T 



2 f 1 1 - (e^) (^ ) + « 1^ (^ ) I 2 ) T ^ 

-2-kS -2tt5 



1/2 



2nS 



1/2 



/ \Un (Z ) ~ ^ n (zo)\ dz 



< 



-2ttS 



< ^Att5 max 1 1 - (s^ (z) + a \U n (z)\ 2 ) I p [U n (z) , tf n (z)] < 

\z\<2ir5 

< ^Aird max [1 1 - e W (z) I + \a\ \U n (z) \ 2 ] p [U n \z) , ^ n (z)} < 

I 2 I <27T<5 



< £p47r<5 max 

ZL \z\<2tt8 L 



l-e^(z)\ + \a\\U mc (z)\^ p[U n (z),* n (z)] . 



(56) 



The last inequality in ([56]) is obtained taking into account the condition max \U n (z)\ < max \U mc (z)\ 

\z\<2k& \z\<2ttS 

which holds for all n — 0, 1, 2, . . . and directly follows from the inequality max \U (z)\ < max \{J mc (z)\ 

J \z\<2tt5 1 wi _ \z\<2ttS 1 v n 

due to fflBT) . We see that, according to (156]) . in the case under study of weakly nonlinear approximation 
([5]) when 



max I" \a\ \U (z)\ 2 ] < max 

\z\<2nS 1 \z\<2nS 



Id \U l 



< max £ 

\z\<2ttS 1 



(i) 



(57) 



the iterations defined by the first equation of (]54]) converge to the unique solution determined by f]M]) if 
the term in the last inequality of (156]) multiplying the norm satisfies the condition 



K a r 

— -47ro max 

2T \z\<2nS 



1 -£ (L) (z)| + \a\ W 



< 1. 



In view of the expression for the transverse wavenumber T = {k 2 — [k sin (</?)] 2 } ^ rewrite the last 
inequality as 



k 2n5 max 

\z\<2tt8 



1 («)| + \a\ \U mc {z)\ <cos(<^). (58) 

Note that according to (!5?1) . condition (|58]) can be written in the form (|37]) with p = as 

g (^-l + «a 2 ) < 1. (59) 



Observe also that (1581) yields the sufficient condition (IM]) for the existence of solution to nonlinear 
IE (125]) obtained in the previous section because the latter reads (E — 1) q < 1. 

We have proved the following statement which constitutes a sufficient condition for the existence of 
solution to nonlinear IE (128]). 



Theorem 6 Assume that the weakly nonlinear approximation (51) holds. Then nonlinear IE ( TI5]) has 
the unique continuous solution if condition j[58\) holds. This solution can be obtained using both the 
iterations defined by the first equation of |5^] ) and the equivalent iteration scheme according to the second 
equation of (5~J\ ) if to consider its solution ^ n (z) as the n+1 approximation (setting \I/ n (z) = U n+ i (z)) 
to the sought for U (z) . 



4 Conclusion 

We have proved, subject to certain sufficient conditions, the unique solvability of the problem of diffrac- 
tion of a plane wave by a transversely inhomogeneous isotropic nonmagnetic linearly polarized dielectric 
layer filled with a Kerr-type nonlinear medium. The diffraction problem has been reduced to a cubic- 
nonlinear IE of the second kind. Based on the use of the contraction principle, sufficient conditions of 
the IE unique solvability have been obtained in the form of simple inequalities. The method presented in 
this work can be generalized so that it will enable one to obtain eigensolutions and soliton-type solutions; 
eigenvalues, also as functions of the problem parameters; and to develop the techniques to wider classes 
of nonlinearities B and operators L of singular semilinear BVPs L (A) u + aB (u; A) = /associated with 
the problems of wave scattering and propagation. 

One the basis of these solution techniques and the IE obtained one can perform numerical investigation 
of the resonance effects caused by certain nonlinear properties of the object under study irradiated by an 
intense electromagnetic field. In particular, one can determine the critical limits of the excitation field 



intensity that govern applicability of the developed mathematical model. The proposed methods and 
results of computations can be further applied to the analysis of various physical phenomena including 
self-influence and interaction of waves; determination of eigenfields, natural (resonance) frequencies of 
nonlinear objects, and dispersion amplitude-phase characteristics of the diffraction fields; description of 
evolution processes in the vicinities of critical points; and to the design and modeling of novel scattering, 
transmitting, and memory devices. 

Appendix 

Let us prove that IE (j251) is equivalent to BVP f|T3|) - f|T6|) ; namely, if U (z) is a solution to IE (1251) then 
E x (y,z) = U (z) exp(i(f)y) is a solution to (fT3l) - (fT5l) subject to representation (}T6|) and vice versa. To 
this end, let us show that IE (125]) and fjl~3]) — fp76]) are reduced to the determination of the solution to one 
and the same BVP and both problems are equivalent to one and the same IE. Indeed, write IE ( 1251) for 
the points \z\ < 2it8 inside the nonlinear layer in the form 

U (z) + y£[F + {z) + F_{z)} = U inc {z) , \z\ < 2ix5 , 
F+(z)= J exp[iT(z-z Q )][l- (eW (z ) + a \U (z )\ 2 )] U (z ) dz , 



-27n5 
2tt8 



(60) 



F_(*) = J expHr(z-zo)] [1- (^ (L) (z ) + a \U (z )\ 2 )] U (z ) dz . 



We have 



U (2tt5) = U inc (2ttS) - i^F + (2tt5) = a inc + a scc; 



2F 

U (-2tt5) = U inc (-2tt5) - ^F„ (-2tt5) = a inc e iinS - l -§rF_ (-2n6) = b scat , 



(61) 



where 



a' 



2 

cat _ _ l J^ F+ ( 27r( J) j b scat = a inc e ^5 _ i^p_ (_ 2vr< 5) (62) 



2 p - -r v n 2T 



denote the quantities (constants) expressed in terms of the solution to IE (1251) . 

Differentiating two times with respect to z the first equality ( 1601) (or IE ( 1251) ) involving functions 
F± (z) and using the continuity condition ()14p for the tangential components of the total diffraction field 
on the permittivity break lines z = 27r5 and z = —2tt5 and representation (ITBI) we arrive at the problem 
in the differential form equivalent to IE (160|) (or (1251) ) 



! r (U)+g(U) = U"{z) + {r 2 -k 2 [I - (e^(z) + a\U(z)\ 2 )]} U (z) = 0, \z\ < 2n5, 
U (2n5) = a scat + a inc , U (-2n5) = b scat , 



(63) 



where the linear differential and nonlinear operators 

£r(U) = U" (z) + T 2 U (z) , g(U) = -k 2 [l - ( e W (z) + a \U {z)\ 2 )] U (z) . 

Note that U (z) in (1631) satisfies the conditions 

U' (2tt5) = iT (a scat - a inc ) , U' (-2irS) = -iTb scat , (64) 

Indeed, differentiating with respect to z the first equality (160|) and setting z = 2ir5 and z = —2tt5 we 
obtain 

U' (2tt5) + i T^F + (2ti5) = -i Ta inc , 

U' (-2tt5) - i (-2tt<J) = -i rV nc exp {AiTn5) , 



which, together with ( I6T1) and ( J62l) . leads to (J64I) . The same result is obtained if we note that on the lines 
z = 2n5 and z = —2tt5, U (z) and its derivative coincide, according to (ITBT) and the continuity condition, 
with the boundary values on these lines of the respective functions U+ (z) = a mc exp { — iT (z — 2 n5) } + 
a scat eX p {iY(z — 2n5)}, U-(z) = b scat exp { — i Y (z + 2 rr 5) } and their derivatives. 

Excluding in fl63|) complex amplitudes a scat and 6 sca * we obtain a semilinear BVP of the Sturm- 
Liouville type 

U" {z) + { T 2 - k 2 [1 - (e^ (z) +a\U (z)\ 2 )] } U (z) = 0, |z| < 2vr<5 , 

iTU (2vr5) - [/' (2tt5) = 2iTa inc , ' (65) 

zr?7 (-27r<&) + U' (-2tt5) = 0. 

We can show independently that (165!) is equivalent to IE (|28l) . Indeed, using Green's function 

G(z,z ) = ^ z -^ 

of the linear differential operator £r{U) we obtain IE (1601) (or ( BHD ) ^y inverting £r{U) in (|63ll with the 
help of Green's function G (z, z ) (that is, by reducing fl65l) to an equivalent IE [5J). The solution to IE 
(j28p satisfies the boundary condition of BVP fl65|) at z = ±27r5 which is verified directly, as above, by 
differentiating with respect to z and setting z = 27r<5and z = —2tt5. 

Assume now that U (z) is a solution to IE (T28]) continuous in the closed interval \z\ < 27r<5 (we note 
that the unique solvability of IE (|28p is proved in Section 4 subject to the sufficient conditions formulated 
in Theorems 1-4). Then constants a scai and b scat are determined from flBTl) and flB^j) so that U (z) satisfies 
boundary conditions in (1631) and (!65|) . and the solution to (fT3|) - (fr5|) is represented in the form ([16]) with 
these constants (which also enter boundary conditions in (|65l) ). 

The same BVP (1651) in the interval \z\ < 2tt5 is obtained from the initial BVP (jT3l) - (TT5l) and rep- 
resentation (|T5|) . This follows directly if we substitute E x (y,z) = U (z) exp (i y) into equation (TT5I) 
taking into account the relationship T 2 = k 2 — cf 2 and the continuity of the tangential components of the 
total diffraction field on the permittivity break lines. 

This statement completes the proof of the fact that IE ff28l) is equivalent to BVP f|T3l) -(fT6l ) . 
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